Quasiclassical Surface of Section Perturbation Theory 
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Perturbation theory, the quasiclassical approximation and the quantum surface of section method 
are combined for the first time. This gives a new solution of the the long standing problem of 
quantizing the resonances generically appearing in classical perturbation theory. Our method is 
restricted to two dimensions. In that case, however, the results are simpler, more explicit and more 
easily expressed visually than the results of earlier techniques. The method involves expanding the 
'phase' of the wavefunction in powers of the square root of the small parameter. It gives explicit 
WKB-like wavefunctions and energies even for certain systems which classically show hard chaos. 



PACS: 05.45.+b, 03.65.Sq, 72.15.Rn 
I. INTRODUCTION AND SUMMARY 

A. Generalities 

The relationship of a quantum system and its classical 
limit has been a major theme in physics for nearly a cen- 
tury. Quasiclassical, or WKB approximations to quan- 
tum theory helped elucidate the relation of the full the- 
ory to the classical limit, but this^was largely restricted 
to the integrable case. Van VleckEJ applied this method 
more generally, but there was little work going much, be- 
yond the integrable case until 1970 when Gutzwillercl in- 
troduced his trace formula, [GTF]. 

Shortly before this, significant advances had been 
made in classical mechanics. Notably, Kolmogo- 
rov, ArjHol'd and Moser proved their famous KAM 
theoremHij, which we discuss further below. Related to 
this and equally important, new ideas and related nu- 
merical studies made chaos an important field of study. 
In practice, the GTF is applied only to 'hard' chaos sys- 
tems, where all classical orbits are unstable. A related 
trace formula, due to Berry and Taboro, applies to the in- 
tegrable case. However, a very large class of systems, the 
'mixed' chaos case, has not been so successfully treated 
by trace formulas. 

It should be appreciated that the trace formula, by 
itself, does not solve the quantum hard chaos problem. 
At best, it gives energy levels, not wave functions. Only 
a few of the lowest levels are even approximately given, 
however. It is mainly convenient for the study of cor- 
relations of widely separated energy levels, which is a 
short time property. However, it is certainly true that 
two decades of efforts to understand the GTF and to 
improve it resulted in techniques which did successfully 
solve the energy level problem at the, WKB level. One of 
those techniques, due to Bogomolnyfl, is the basis of the 
present work. This niethod is closely related to that of 
Doron and SmilanakyQ, and more subtly to the ideas, of 
Berry and KeatingQ and of Cvitanovic and EckhardtQ. 



B. Perturbation theory 

It is not necessary to emphasize that approximations 
based on large or small quantities are important in math- 
ematics. Indeed, KAM theory is a resolution of prob- 



lems of classical perturbation theory first uncovered by 
Poincare. Attempts have been made, with some success, 
to quantize a generically perturbed classically integrable 
system. Of course, until KAM the classical system itself 
was poorly understood. However, KAM basically found 
that a perturbed system generically looks rather like a 
mixed chaos system, and we have remarked that quan- 
tum mixed chaos is far from being understood. 

On the other hand, given an integrable quantum sys- 
tem, perturbation theory on it is routinely carried out. 
This is relatively easy for the case when h is not 'small', 
that is, when the perturbation is smaller than the energy 
level spacing. 



C. Two small parameters 

In this paper we are interested in the case where two 
small parameters are present. The first parameter, de- 
noted e, gives the scale of the difference between an ex- 
actly solvable, integrable case, and the 'perturbed' case 
of interest. The second parameter is the dimensionless 
Planck's constant, h, which gives the scale of the leading 
order quantum or wave effects. 

The relationship between these parameters is crucial, 
and the theory is exceeding rich as a result. Although this 
has certainly been understood for some time, we have not 
found a very clear discussion in the context of the formal 
perturbation-quasiclassical approaches used up to now. 

The relationship comes from the well known fact, made 
manifest by Feynman's path integral formulation, that 
quantum effects smear classical phase space structures 
over areas of size h = 2-Kh. We choose to compare 
Planck's constant with the area of the leading phase space 
structure, call it S. The dimensionless Planck's constant 
is h/S, which we continue to call h, i.e. we choose iS 
as the unit of action. For example, in a nearly circu- 
lar billiard of radius i?, S could he S = pR = hkR, 
so h/S = 211 /kR = X/R. Here p is the momentum, k 
the wavenumber and A the wavelength of the particle in 
the billiard. If the dimensionless Planck is of order unity, 
only the gross features of the classical system are reflected 
in the quantum properties. But for small e, or in the pres- 
ence of chaos, there are typically classical structures on 
all scales, i.e. on the scale of e^, where /3 ranges from 
zero to infinity. Typically, the smallest nonvanishing /3 
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of interest is /? = ^. If /h is of order unity, it means 
that quantum systems do not reflect the structure at the 
levels e'^ with [3 > M but are sensitive to structure with 
P<M. 



actions are regarded as functions oi 9 — 9' = 2TruJi/uj2 
and E. We note that S is just the integral J pdq along 
the trajectory from one surface of section crossing to the 
next. 



D. Integrable classical systems 

An integrable classical system has enough constants 
of the motion that each classical orbit in 2c? dimensional 
phase space lies on a d-dimensional surface in that space. 
These surfaces turn out to be tori. The d constants 
of the motion can be taken to be the action variables, 
Ij, j = 1..0?, and they label the invariant tori. From 
now on, we take d = 2. In these canonical variables, 
the Hamiltonian Hq depends only on the actions, i.e. 
Hq{Ii,I2) ~ E, where E is the energy. The conju- 
gate angle variables, ^1,^2, satisfy Hamilton's equations 
9j = dHo/dlj = ojj = ujj{Ii,l2)- The angle variables 
fix the point on a given invariant torus. We shall also 
assume, to avoid certain complications, a principle of 
'sufficient nonlinearity', which posits that the winding 
number lui/uj2 is not too constant as Ii and I2 are varied 
at fixed E. We also often use two other constants equiv- 
alent to the /'s, usually E and some other convenient 
variable, such as the winding number. 

A common way of describing the perturbation is to 
give a Hamiltonian 



H ^Ho{hj2) + eHi{h,l2,0i,92) 



(1) 



where the perturbation is periodic in the angles. It is well 
known in this field, however, that other formulations are 
often more convenient. 



E. Poincare surface of section 

A Poincare surface of section, SS, is a phase space 
surface through which all tori of interest pass, once and 
only once. In action angle variables, or their equivalent, 
we usually take the surface 92 — 0, E fixed. Fig. 1 
illustrates such a surface. A given orbit crosses the SS 
each time 6*2 is a multiple of 27r. Given 9'^ and I[ when 
92 = 0, the dynamics predicts a new (/, 9) = T{I' , 9') at 
the next intersection of the orbit with the SS. (From now 
on, we call the variables on the SS simply 1,9, without 
subscripts.) T is called the surface of section map and of 
course it depends on E. For an integrable system, 1 = 1', 
9 = 9' + 2-ku)i/lo2- For a non integrable system, / is not 
constant. 

The surface of section is usually displayed in a graph 
of / vs 9, for Q < 9 < 2tt, and with points 0, 27r identified. 
The intersection of the invariant torus / = const with the 
surface of section, is the horizontal line / = const on this 
graph. 

The surface of section map is conveniently given by a 
generating function. This is an action, S{9,9') such that 
/ = dS{9,9')/d9, r = -dS{9,9')/d9'. In the integrable 
case, S{9, 9') = Sq{9 - 9') = {9- 9')h + 27r/2, where the 
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FIG. 1. Schematic of the definition of surface of section 
62 = const. An orbit starts at 6\ = 9' in the SS, goes around 
the hole once, arriving aX 6\ = 6. It advances by an angle 
6-6' = 2-KUJi/uj2 on the SS. 

Under perturbation, S becomes 

S{9, 9') = So{9 - 9') + tSi{9,9') + 0{e^). (2) 

There is no difficulty, except tedium, in using perturba- 
tion methods to find Si given say. Hi, since the calcula- 
tion involves only short orbits. That is to say, the funda- 
mental difficulties with classical perturbation theory lie 
in understanding the long time behavior. 

F. Quantum surface of section method 

Our theory is expressed in terms of S. Bogomolny^ 
has shown (without requiring e small) that the surface of 
section transfer operator or kernel 



T{9,9') = 



/ 1 d'^S{9,9') 
\ 2mn d9d9' 



exp( -5(0,0') ) (3) 



can be used to find the energy levels of a system, in quasi- 
classical approximation, although generally the solution 
is not very explicit. We mention some of the main results 
below. Bogomolny's method is by no means restricted to 
action angle variables or ajw, particular surface of section. 
There is a large literatureil't3 exploiting and verifying this 
technique. We call this the quantum surface of section 
method. We shall find that using Eq. (||) in Eq. (^) and ex- 
ploiting the small parameter e, allows a rather complete 
and explicit solution of the quantum problem. 



G. KAM theory 

We summarize KAM theory as follows. A Hamilto- 
nian H = Hq + eHi is assumed, where Hi is a suitably 
smooth and bounded function of Ij , 9j , and is periodic 
in the angle variables. It is found that, generically, the 
original invariant tori which are rational, or equivalently 
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resonant, are destroyed, as well as those in their immedi- 
ate vicinity. A rational torus is one for which the wind- 
ing number uji/lo2 is a rational, p/q. Any orbit on such a 
torus whose 62 increases by 27rg will have 9i increase by 
2-Kp : itjis thus a periodic orbit. The Poincarc-Birkhoff 
theoremU says that only a few, often just two, such or- 
bits of the original continuum on the torus, survive the 
perturbation. One surviving orbit is stable, the other 
unstable. 

Within a width in action, \fel-pq about the rational 
torus, the neighboring tori, both rational and irrational, 
are destroyed or modified strongly. Here /pq is a char- 
acteristic classical scale of action associated with the pq 
torus. 

The rational tori are of course dense. However, the 
series I^q converges to a finite action, if the pertur- 
bation is smooth enough, so the total volume in which 
the original invariant tori are destroyed is a small but 
finite fraction of phase space, proportional to ^JT^. 

We may 'cut out' this phase space volume near the 
original rational tori. What is left is called a KAM sem. 
The KAM set is a Cantor set, which however occupies 
most of the volume of phase space. The tori in the KAM 
set are invariant and only slightly modified under pertur- 
bation. This is the main result of KAM theory. 

The destroyed tori are largely replaced by new invari- 
ant tori. These have a new topology, namely the topology 
of tubes with the surviving stable periodic orbit (s) at the 
center(s). In the surface of section representation of Fig. 
2, these look like islands. 




ratrix, there are invariant tori, perturbed on the scale of 
^/elpq which 'gradually' go into the KAM set. 

Fig. 2 shows such intersections of invariant tori with 
the surface of section, which we call invariant loops. 
These represent tori in the KAM set, new invariant tori 
and separatrices. 

The new invariant tori can also be classified as rational 
or irrational, and the whole process must be repeated, ad 
infinitum, giving rise to secondary resonances, and so on. 
However, the original perturbation scale e, is replaced for 
the secondary resonancesa by something like e^/"^, which 
is smaller than any power of e. p. 

The scale of the stochastic region is also very smalla, 
namely of order exp(— l/-ye), again beyond any ordinary 
perturbation theory. These phenomena are represented 
schematically in Fig. 2. 

KAM theory thus shows that after perturbation phase 
space consists of an infinite number of parts. 

1. A large part, the KAM set, with a Cantor set struc- 
ture, of relative phase space volume of order 1 — ^/e 
but weakly and smoothly perturbed on the scale of 
e. 

2. A relatively big part, consisting of new invariant 
tori near the original rational tori, which usually 
has a scale of ^/e, but made up of an infinite num- 
ber of pieces of size y/elpq. Remark that ^/e is big- 
ger than any power of e, which is one reason that 
perturbation theory in e must fail. 

3. Conceptually important but smaller-than-any- 
power-of-e stochastic regions which also are near 
the original resonant tori. 

4. The still smaller secondary resonances and stochas- 
tic regions, near the new rational tori. 

5. The smaller yet tertiary resonances and stochastic 
regions, ad infinitum. 

A final remark for this section. The appearance of 
the square root -^/e" is generic in the sense that it occurs 
when the second derivative Sq{9) does not vanish at a 
value of Qpq = 2'Kp/ q explained further below. It is easy 
to construct examples for which this fails. Our method 
can be generalized to deal with such cases. 



FIG. 2. Surface of section invariant tori and separatrices, 
showing the scale \/t of resonances, the scale e of the 'golden 
mean' torus in the KAM set, and secondary resonances, all 
for e = 0.01, r? = 0.19. 

Separating these tubes from one another and from the 
KAM set is a separatrix region, consisting of the stable 
and unstable manifolds of the remaining unstable peri- 
odic orbit. Because these manifolds intersect transversely 
at a homoclinic point, although at a very small angle, this 
separatrix region becomes a chaotic, homoclinic tangle, 
also referred to as a stochastic region. Outside the sepa- 



H. Results 

Most of these results have been announcedS. We are 
able to find explicitly all energy levels and 'wave func- 
tions' to leading order in fi. By wave functions, we mean 
eigenstates of the T operator, which arc closely related 
to the usual wave functions. The equation we solve is 

TV- = V' (4) 

This equation can be solved only if E, on which T de- 
pends parametrically, is on the spectrum or more pre- 
cisely is on the WKB approximation to the spectrum. 
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There generally exists an exact operator or kernel K to 
which T is a quasiclassical appMximation, for which this 
procedure gives exact answeraly. There is a Fredholm 
integral equation over the SS 

r) = u{9', r)+ j d9K{9', 9; E)^i{9, r). (5) 

and the Green's function G{r',r;E) can be computed 
from the solution /i. Eq.(U) is the homogeneous version of 
Eq.(||) with K approximated by T. For billiards, with the 
boundary as SS, Eq.(|^) is the equation for the boundary 
integral method, and ip is the normal derivative of the 
usual wavefunction ^'(r) on the boundary. 

There is a known procedure to find ^' given ip. How- 
ever, contains most of the information desired rather 
directly, and it is not necessary as a rule to find In 
fact, there are many ways available to represent a com- 
plex function of two variables, and -0 is a useful one even 
in nonperturbative contexts. 

We can describe what we have done as the solution of 
Eq. in WKB approximation for the wave functions and 
energy levels as a power series in y^, or in e, depending 
on the phase space region. We also have a criterion for 
whether the series should be ^/e, or e. Thus, in principle, 
we find explicit solutions for those quantities which are 
not smaller than any power of e. To be a bit more ac- 
curate, we find explicit solutions for the leading ^ and 
first few powers, and provide a procedure which rapidly 
becomes tedious, for the higher orders. In this paper we 
will concentrate on the leading order or two, and not on 
issues of the breakdown of the series, etc. 

A second interpretation is in terms of quantum per- 
turbation theory. We do a quantum perturbation the- 
ory about an initial state ip = exp(iZ0). However, unless 
y/c/h << 1, this is not an ordinary perturbation theory, 
but rather a degenerate perturbation theory, in which of 
order y/e/h unperturbed states are mixed, in the reso- 
nant case, 01 e/h states in the KAM set case. This can 
be seen in Fig. 2, when it is realized that unperturbed 
states are represented by action values spaced h apart. 
So, we can effectively carry out the diagonalization of the 
matrix mixing these states in this case. 

The fact that a form of degenerate perturbation theaty 
is needed was apparently first pointed out by Shuryakllj, 
in the case of resonances. Thus, if \/eIpq << h ordi- 
nary quantum perturbation theory works well. Above 
the 'Shuryak border', ^/elpq > h, ordinary perturbation 
theory breaks down as a number y/elpq/Ti unperturbed 
quantum states are strongly mixed by the perturbation. 
However, even in the case of the KAM set, if e/?l is large, 
many unperturbed levels must be mixed to obtain a so- 
lution. 

A third interpretation of what we do is that we are able 
to find explicit approximate formulas for the invariant 
loops. These loops have the topology of circles and are 
invariant under the surface of section map. The loops 
are found as a power series in e for tori in the KAM set. 
We find explicit formulas for the new, resonant, invariant 
loops as a power series in ^/e. We find explicit formulas 
for the separatrices in this way also. 



The loops can be quantized, in essence by EKB theory, 
and the formulation then gives the other quantization of 
the motion away from the SS. 

We also have expressions for the scale of the hole 
around a resonant torus, Ipq. Having these formulas, we 
can find the quasiclassical quantization conditions and 
the wave functions. 

The series in e or ^/e' that we find are not expected 
to converge, but are probably asymptotic. If they con- 
verged, we would have found exact new constants of the 
motion, which cannot exist. 

We, however, do not find any stochasticity, or any 
secondary resonance phenomena, since these are smaller 
than any power of e. This has the unfortunate impli- 
cation that our work does not shed much light on the 
quantization of mixed chaos. It would indeed be surpris- 
ing, however, if the quantized wave functions of extended 
stochastic regions had simple analytic expressions. In 
fact, mixed chaos consists of three regions, a KAM like 
region with invariant tori, a region presumably strongly 
chaotic and treatable by the GTF and its improvements, 
and a transition region between the two. It is the tran- 
sition region which has resisted attempts to quantize it. 
There is some workE3 bearing on this transition region 
which exploits the fractal nature of the region. 

Our results are simple and familiar on the one hand, 
but rather complex on the other. That is, we obtain re- 
sults of a type familiar from standard WKB theory. But, 
the results are complex in the sense that the results of 
classical KAM theory are complex. Various cases are 
possible depending on the choice of system. Since the 
quantum theory must follow the classical up to a certain 
level of approximation, depending on h, the quantum re- 
sults must also be complex in this sense. 

I. Other work on combined quasiclassical and 
perturbation theory 

Combined QCA and PT has been used and studied 
extensively. The general scheme is to find approxima- 
tions for the invariant tori and then use EKB methods 
to quantize them. In this sense our approach is a varia- 
tion of this already existing work. 

The difference is that we need to find only the invariant 
loops in the surface of section, and not the whole torus to 
which an invariant loop corresponds. This has the usual 
advantages of using the surface of section. It reduces 
the dimension, thus simplifying calculations and making 
it easy to visualize the results. Our results are rather 
simple to derive and write down. The older approach, 
in practice makes approximations for the invariant tori 
which are completely numerical or involves expressions 
with many terms requiring computer algebra. 

An important techiuque is that of the Birkhoff- 
Gustavson normal fornJlJlil. This is closely related to 
the Darling-Dennison formll3. In this case, the starting 
point is two or more harmonic oscillator Hamiltonians 
coupled by terms consisting of polynomials in the dis- 
placements and momenta. Successive canonical trans- 
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formations eliminate the coupling terms, order by order, 
in favor of a normal form Hamiltonian, which basically 
consists of products of powers-of the original harmonic 
oscillators. Ozorio de AlmcidaHl has used this to obtain 
results somewhat similar to ours. 

This method is particularly relevant in treating vi- 
brational modes of small molecules. Our method needs 
modification for this case, as the unperturbed winding 
number is constant. On the other hand, it's not clear 
whether Birkhoff-Gustavson can be generalized to deal 
with weakly distorted billiards, for example. 

In practice, Birkhoff-Gustavson relies on computer al- 
gebra and numerical evaluations. Although discussion of 
the convergence of the series is givenllj, mention is not 
made of the relationship of the small parameter to h. 

Another approach is to propagate numerically an orbit 
which lies on an invariant torusHS. Fourier methods then 
allow a determination of the fundamental frequencies, the 
Wi's and then the torus itself. This method is basically 
numerical, but it is quite efficient. 

Although the SS technique is often used to display 
results, it is seldom used as a calculational tool. The 
one exception we have found is the work of Noid and 
Marcualj. Indeed, they used two surfaces of section, SSI 
and SS2 which cut the invariant tori in topologically dif- 
ferent ways. By numerical propagation of an orbit, they 
found the invariant loops li and I2 on these two SS. EBK 
quantization makes the area of the loops take the form 
areai = h{ni + i), area2 = h{n2 + assuming the 
usual i for librational states. In two dimensions, the SS, 
and thus in effect the orbits, depend on two parameters, 
one of which is taken as the energy. These parameters 
are varied until the areas of the loops have the above 
relation to integers. Noid and Marcus obtained numeri- 
cally good results for the energy levels. Use of SS's gave 
some advantage computationally. The method does not 
rely on a small parameter, but only on the existence of 
sufficiently simple invariant 



A different kind of effort 



studies the transition 



away from the Berry- Tabor trace formula, valid for com- 
pletely integrable systems at e = 0, and towards the 
Gutzwiller trace formula, usually considered only for 
hard chaos. The Shuryak border and the are not 
mentioned in this work, and an expansion is made in e. 
This seems peculiar, at first sight, since the trace formu- 
las are purported to give the density of states. 

In fact, to get energy levels from the trace formula, 
a resummation of the divergent series must be carried 
out. In other words, the contribution of very long orbits 
must be somehow included. One method for doing this 
uses Bogomolny's T and is basically equivalent to Eq.(||). 
The trace formula without resummation is directly useful 
if an appropriate energy average is made, as it predicts 
energy-level correlations over a correlation distance long 
compared with the mean level spacing. This energy aver- 
age is equivalent to looking at short time orbits, for which 
ordinary perturbation expansions in powers of e are ad- 
equate. The work based on the trace formula restricts 
considerations to this case. We can readily reproduce the 
cited results and the trace formulas above by use of the 



transfer operator. 

Tiere are also rigorous mathematical results in this 
fielcD. We are not sure how to make a simple statement of 
what has been rigorously proved. It seems that Lazutkin 
has succeeded in quantizing the KAM set at the EKB 
level, and has shown that the number of levels found is 
correct. The wave functions are also correct, up to the 
possibility discussed below of accidental degeneracy. 



J. Applications 



There are a number of papera2Jc3 in the recent litera- 
ture which have uncovered interesting phenomena in this 
perturbative case. Indeed, it was this literature which led 
us to develop this theory. Some of this work has prac- 
tical applications in the construction of tiny lasers with 
almost circular resonant cavitiescil. 

The other work is more concerned with localization 
theory. The existence of a rather complete and analytic 
solution was not suspected, so statistical studies of local- 
ization and energy level statistics were undertaken with 
supporting heuristic arguments. 

In these papers, the basic object is a weakly deformed 
circular billiard, which can be described by a boundary 
in polar coordinates R{d) = Rq + eAR{6). In these bil- 
liards, the role of h is played by fc~^, the inverse dimen- 
sionless wavenumber or square root of the energy. Fig. 
3 shows such a billiard together with some characteris- 
tic lengths assoeiated with the perturbation theory. The 
older literatureESl treats such boundary perturbations in 
more ordinary quantum perturbation theory, valid only 
for k^/e « 1. j-jp 

The case AR = cos some special features in 

that a nearly circular elliptical billiard has this as its lead- 
ing order term, and such a billiard is integrable. Higher 
order corrections are quite interesting since they are re- 
sponsible for island chains, etc. That is, this case can be 
thought of as an integrable elliptical billiard with lead- 
ing perturbation proportional to e^. We reserve a more 
detailed analysis to future publication, however. 



K. Beyond KAM theory 

KAM theory comes in several flavors, depending on 
the smoothness of the perturbation. Examples where 
KAM does not apply have also been of recent intprpst 

A nearly circular Bunimovich stadium billiardcjO'EHl 
has AR ~ |sin6'| — 2/tt. This gives the leading order in 
e deviation from a circle. We call it the nearly circu- 
lar stadium billiard. The constant 2/7r is subtracted to 
make the average perturbation vanish, or equivalently, 
to make the perturbed billiard have the same area as the 
unperturbed circle. We shall usually drop it for brevity, 
however. For comparison purposes we shall generalize it 
to AR — \/ sin^ 6 + ?f to give a billiard with a similar 
shape, but with a smoothing parameter rj. We call this 
the 'nearly stadimn' billiard. 
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It is rigorously knownEa, for 77 = 0, that because of 
the lack of second derivative, the long time behavior 
is diffusive in angular momentum space, and there are 
no invariant tori at all. If rj is finite but small, say 
•q « -y/e, similar diffusive results are obtained numeri- 
cally. A somewhat related case is to take Ai? as a 'ran- 
dom' finite Eaurier series constrained to have M maxima 
and minimao. There is a third, large, parameter M or 
1/77, whose relationship with k and e must be specified. 
The possibility of large higher derivatives of the pertur- 
bation can modify or invalidate KAM theory as well as 
stationary phase methods. 

The general interest in the quantum version of these 
non KAM cases has been to find changes of wave func- 
tion localization behavior as a function of the parameters, 
and to interpret these regimes as belonging to different 
types of dynamical localization. At some large enough fc, 
the large derivatives in these cases start to play a role, 
and stationary phase arguments and perturbation theory 
breaks down. The contribution of our paper is for values 
of k less than this characteristic mamitude. A detailed 
study of this breakdown is underwayEll. 

However, to some order in -y/e, for the resonances, ap- 
proximate invariant tori continue to exist, and we can 
find them and quantize them. If these tori break down 
at order r, then for ke^^"^ « 1, the quantum system 
acts as if the tori exist. For larger fc, quantum mechanics 
can 'see' the dissolution of the tori, and our theory no 
longer applies. This transitian .ruid its interpretation is 
the theme of several papersOE^Ej'EJ. 

This can also be interpreted as the fact that quan- 
tum mechanics and its limiting classical mechanics may 
have quite different long time behavior. The two approx- 
imately coincide up to a time t^, a time usually inter- 
preted as ?i/ A where A is the mean level spacing of the 
quantum levels which are appreciable in the initial wave 
packet. Beyond this time, the classical behavior is ir- 
relevant to the understanding of the quantum behavior. 
Thus, the long time classical behavior, which is the fo- 
cus of most classical theory, can be irrelevant unless it 
also describes the short time behavior, t < Tn. Quantum 
mechanics does not 'care' if invariant tori are indeed in- 
variant. It only cares that the classical orbit stays close 
enough to a torus for long enough that phase interference 
effects can establish the quantum states— 

A buzz word used in this contextESEJ is 'cantorus', 
which is a Cantor set invariant phase space structure re- 
placing an invariant KAM torus whids. has disappeared 
as a parameter is changed. CantoriE3 can 'trap' for a 
while classical orbits, thus influencing the short time be- 
havior. Although this concept is well established in cer- 
tain contextsL^I, it is not clear to us exactly where the 
cantori arc and where they are not in the present case. 
In any case the relevant short time behavior can be de- 
scribed much more simply and cogently than by invoking 
cantori. The invariance of a cantorus is what makes it 
interesting and difficult in the classical context, but this 
invariance is of no relevance for quantum systems. 

There are also systems which have no explicit small 
parameter. Nevertheless, there can be classes of eigen- 



states which can be obtained by our technique. The 
best known of these are the whispering gallery modes 
in smooth convex billiards and the, quantized 'bouncing 
ball' orbits in the stadium pbilliarcO. There are also spe- 
cial modes found in workc3 on a ray splitting billiard, 
and this was the first application of the methods of this 
paper. To support a quantum state that we can describe 
in our theory, there have to be many similar orbits, e.g. 
nearly nonisolated orbits similar to the orbits on a KAM 
torus. Basically, a small parameter appears as the solu- 
tion is found. The potentials found below (which are of 
order ke) are in this case replaced by a big potential, but 
only 'low energy' solutions are accurately found. 



II. OUTLINE OF THE METHOD 
A. Specific model system: Perturbed circle billiard 

For a perturbed integrable system, it is possible to 
choose coordinates and a surface of section so that the T 
operator, Eq.(^), is based on an action S, of Eq.(^). To 
be explicit, we consider a perturbed circular billiard. 

In quantum language, we take units Rq = 1, h — 1, 
particle mass = 1/2, so fc = ^/E is the dimensionless 
wavenumber, equivalent to l/h. We take the billiard 
boundary OB as SS. Then Bogomolny's unitary opera- 
tor isQ 



T{e,e',k) 



/ k d^L{e,e') 

\ 2m dedO' 



exp {ikL {0, e')) (6) 



where L is the chord distance between two points, spec- 
ified by polar angles, on OB. Expanding, 



kL{e,e') = 2k 



sm ■ 



1 + e 
•■) 



Ai?(6i) + AR{e') 



(7) 



We assume that k is large, and e is small. We introduce 
6 — \ft. The relation between k and h is specified by an 



» 



integer M, such that nominally kb ^ 1, and kb 
1 while kb'^'^^^ « 1. In phases that are expanded in 
powers of b, we must keep terms up to kb^^ , even though 
they are much smaller than other contributions to the 
phase. This is because such phases give changes of order 
unity to the wavefunction as the angle is varied. In this 
paper we shall usually consider M = 2, i.e., we shall 
consider kb^ << 1. 

In prefactors, we may keep only the leading order 
terms and succeed in making only fractionally small er- 
rors. In finding the prefactor of the wavefunction we shall 
need to keep terms of order unity, which we may regard 
as being of the same order as 

This enables us to make a convenient, but not strictly 
necessary, approximation to L2 above. We will later find 
that only values of 9— 9', near a fixed angle 6/, depending 
on the angular momentum range considered, are impor- 
tant. In L2, we replace sini(6' — 9') by the constant 
d = |sin ^6; I , an approximation accurate if kb^ is small. 
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B. Bogomolny Integral Equation 



C. Simplest case 



According to Bogomolny, the energy levels of the sys- 
tem are given in QCAlI by solutions E = Ea = of 

D{E)^dct{l-T{E)) = 0. (8) 
This equation is usually approached in one of three ways. 

1. The imaginary part of the logarithmic derivative of 
D, dln{D) / dE can be expanded in traces of powers 
of T, which yields, for example, the Gutzwiller trace 
formula, if all periodic orbits are isolated and unsta- 
ble. If T is integrable, this gives the Berry- Tabor 
result, and in the present case, it gives thcr-aec^ 
turbed Berry- Tabor results mentioned aboveE^TcJ 
The original derivation of these results is more dif- 
ficult, and it had the defect that the sum was math- 
ematically rather ill defined, since it was not abso- 
lutely convergent. The organization of the series 
by the T operator derivation groups together or- 
bits coming from the same power of T and at least 
yields a series which either converges or diverges. 

2. The Fredholm determinant Djiasf/ be expanded by 
the rules of Fredholm theoryQ'E£l, giving an abso- 
lutely convergent expression, which in quasiclassi- 
cal approximation is a finite sum. An important 
improvement uses the unitarity of T to make each 
term in this sum real. The same traces of powers of 
T and periodic orbits appear as in the trace formu- 
las, but organized into 'composite orbits' 
the main result of 'resummation' of the GTI 

3. The kernel T may be represented by a discrete ma- 
trix, and numerically diagonalized. Although this 
method is disdained by those committed to periodic 
orbit theory, it gives very good results. 

Our seemingly more difficult technique finds wavefunc- 
tions ip satisfying 



(9) 



which can be done only for k = ka- Our method is 
tractable only if the orbits are nearly nonisolated, as is 
the case for a perturbed integral system. 
This problem is naturally generalized to 



(10) 



which can be solved for all fc, and k is regarded as a pa- 
rameter. This allows a study of perturbed quantum maps 
in addition to two dimensional autonomous systems. The 
phase Lo will be a function of k and usually, it is rather 
easy to solve the equation u}{k) = 27rn which gives the k 
values such that a solution of Eq.(p|) exists. 



For pedagogical reasons we start with a case special in 
two respects. First, we assume that the states are entirely 
made up of 'low' angular momenta, which means that ij: 
is slowly varying compared with the variation of the T 
operator. From a classical perspective, the relevant orbits 
will pass fairly close to the center of the circle. Straight- 
line orbits of this class, for the billiard of Fig. 3, pass 
through a shaded region at the billiard center. 




FIG. 3. A nearly circular, nearly stadium billiard 
R{e) = Ro + AR{e) with e = 0.3. The straight line seg- 
ments of a stadium have a length of 2e. Its area is equal to 
the area of the dashed circle _Ro = 1- If approximated to the 
first order in e, AR{9) has cusps at 9 = 0, tt. Orbits that 
cross the shaded area are affected by Z = resonance. Several 
relevant lengths are shown. 

Second, we assume that ke << 1. This means that the 
wavelength is long compared to the shift of the boundary 
eAR. Contrary to intuition, as we mentioned earlier, this 
condition does not suffice to make perturbation results 
trivial. For that, it is necessary that ky/e << 1, 

With this condition, we may expand the Eq.(^ as 



T{e,e') 



-L(l + ...)eM2ik-'-^S0^ + ...)x (11) 



X [1 + ikeV{e) + 



(12) 



where S9 = 9' — 9 — tt is smaU, i.e. 6' ~ 9 + tt. [We 
arrange things, to simplify the notation, so that 9' ~ 9+tt 
rather than 9' = 9 - tt.] Here V{9) = L2{9,9 + tt) = 
AR{9) + AR{9 + tt). We then expand 



^Ij{9') = V(^ + tt) + 59^/(9 + tt) + -S9^^j"{9 + tt) + 



We also expand 



(13) 



(14) 



where Em is a shift of eigenphase to be determined and 
luq is defined below. The ^tt comes from the prefactor 
and Gaussian integral over S9. Using expressions (11) , 
(|l|) and ([l|) in Eq.(|l^), and doing the 69 integral, we 
find the conditions for a solution. 
First, we must require 
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(15) 



In this case, luq is or tt, since must be 27r periodic. 
Insisting that the leading small terms vanish, we find that 
Tp must also satisfy the equation 



1 



(16) 



This is a familiar equation, similar to the equation of a 
quantum pendulum, for motion in a periodic potential of 
strength k^e — (fcy^)^ , and unit Planck. Alternatively, 
we may take a potential whose scale is unity and think 
of Planck's constant as h — l/{ky/e). Thus b — -^/e natu- 
rally appears. This is such a well known and thoroughly 
analyzed equation, that the problem can be considered 
solved. 

If kb is sufficiently small, we can make an ordinary 
quantum perturbation expansion about the starting state 
of zero angular momentum ip = constant. If kb is larger, 
some other method must be used. 

If kb is large, standard WKB theory makes the Ansatz 
^{9) = exp{ikbfi9) + ig{9) + 0{l/kb)) where df/d9 ^ 
/' ~ 1 and g ~ 1. Using this approach gives 



f'{9) = ±^Ern - V{9) 



(17) 



which is TT periodic. If Em > V for all 9, we may quantize 
Em by the condition 



kb{f{9 + Tl)~ f{9))=±kb d9y/Em - V{9) 

JO 

= Wo + 27rm. (18) 

We call this the 'rotational' case by analogy with the 
pendulum. 

The usual next order analysis gives g{9) = 
i^\n f'{9), which is customarily written as a prefactor 

{Em-V{9))-'/\ 

If Em < nrd.xV(9), there will be potential wells, at 
least two in this case, with corresponding librational mo- 
tion. This may be treated at various levels of approxima- 
tion. If tunnelling between the wells is neglected, there 
will be a quantization condition 



kb 



d9y/Em - Vi9) = 7r(m + ly). 



(19) 



Here, the limits are the angles where the square root 
vanishes, the classical turning points, where the leading 
WKB approximation breaks down. 

The Maslov index i/, usually i, is included. It can be 
found, for example, by the usual device of approximating 
V hy a, linear function in the neighborhood of the turning 
points, finding the Airy function solution, and using it to 
interpolate between the WKB solutions away from the 
turning point. 

The wavefunction, sufficiently inside the turning 
points, is approximately given by 



4>^iEm~ V) 



-1/4. 



kb 



d9\ E, 



V{9) + I 



(20) 



This gives a two-fold degeneracy, since there are libra- 
tional levels in two identical wells. If desired, the expo- 
nentially small splitting of these levels can be estimated. 

The condition Em = maxy(0) gives the separatrix 
of the motion between librational and rotational motion. 
Again, the simplest WKB approximation in the neighbor- 
hood of the separatrix can be corrected by well known if 
tedious methods. 

Is it also straightforward to take into account sym- 
metries which may exist, for example under 9 < — > —0, 
or time reversal, which guarantees that eigenfunctions of 
Eq.(p|) can be taken real. 



D. Small angular momentum, fee > 1, ke^^'^ << 1 

If fee > 1 it is not possible to expand the exponential 
representing the perturbation in the T operator. We can 
show, however, that the WKB solution of the previous 
section remains valid. 

Expanding the phase of the T operator about the 
point 9' — 9 + TT, as before, we obtain 



S{9,9') ^kL^2k 



^k69^ 



ke{AR{9) + AR{9')) (21) 



[In Eq.(21) we replaced Lq by 2, its stationary value, 
when multiplied by e.] 

The WKB Ansatz is ip = exTpikbf{9). The same pref- 
actor will also be found, but we ignore it for now. 

Returning to Eq.(^), and using this Ansatz, we ex- 
pand all functions of 9' about 9 + tt. I.e. f{9') ~ 
f{9 + tt) + 69 f {9 + tt) to order 69, [since kb « k] and 
AR{9') w A R( 9 + tt), since kb » fee. Doing the integral 
reduces Eq.(|3J) to 

exp[ikbf{9) + icu] 
= i exp [i {2k + [kbff/k + k(y{9) + kbf{9 + tt))] (22) 

where V{9) = AR{9) + AR{9 + tt) as before. 

For Eq.(p2|) to hold, the exponents of order kb must 
combine to give a constant uj^, i.e. f{9 + tt) — f (9) + 
bj^/kb. At order kb^ solution is possible provided (/')^ -I- 
V{9) is a constant, which again we call Em- Thus again 



f{9) - ± / d9'^Em - V{9'). 



(23) 



The lower limit can be chosen at our convenience. Notice 
V{9) = V{9 + tt)=^ f{9 + tt) = f{9)+ const. 

Everything goes through exactly as in the previous 
section for the rotational states where Em > uiaxV. For 
Em < max V we have librational (or near separatrix) 
states where leading order WKB fails in the usual way. 

However, we may proceed by almost the standard 
technique. Namely, for angles near the turning point, 
ke{Em — V) is small, even if ke is not. In this region, 
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therefore, we can expand as in Eqs.(|ll|), (^3|). The solu- 
tion in this region interpolates between the regions where 
WKB is good and we obtain the usual results, including 
prefactors and Maslov indices. 

In Fig. 4 we compare the results of this approximation 
with exact numerical determination of the wavefunctions. 
The exact results show the tunnelling tails and Airy- 
function-like interpolation at the turning points, which 
we have not bothered to calculate in WKB theory. 




71/8 71/4 371/8 71/2 

angle 



FIG. 4. 'Stadium' potential |sinS| vs angle. States and 
potential are symmetric about zero angle. Bound and contin- 
uum, WKB and exact states are shown, with zeroes at WKB 
'energy' parameter Em- k^/i = 42.3 is fixed. Inset: Angular 
momentum representation of continuum state m = 48 and 
exact state near angular momentum m = 168. 

We remark that in this case, the results are indepen- 
dent of the size of fee, that is, there are no corrections 
to the result of order fee, but only of order ke^^^. Thus 
we show two wavefunctions that cannot be distinguished, 
with the same value of kb, but quite different values of 
kb^, one smaller than unity, the other larger. 

We remark that the 69 integral is effectively over a 
width of order I/VX. However, the center of the effective 
window of integration is not at 60 = 0, if we take into 
account the existence of /, but rather at 69 ~ 2y/ef'. This 
shift will be small compared with the width if ^/e << 
yJVjk, or /ce << 1. Thus the term in of Eq.(^ docs 
not contribute for kt << 1. In the present section, that 
condition does not hold, but in the turning point region 
where we are expanding, the shift is small because /' is 
small. 

Thus, we conclude that, even if fee > 1, the solution is 
given by the solution of Eq.([l6|), using the usual WKB 
methods. 



E. Classical interpretation 

We now give the interpretation of these results in clas- 
sical terms. 



The action S{9, 6'), the phase of the T operator, gen- 
erates the surface of section map T{1',9') {1,9) by 
the equations / = dS/d9, I' = -dS/d9'. For the dis- 
torted circle billiard, this is explicitly I = —k{9 — 9' + 
tt)/2 + k€AR'i9); I' = -fc(6' - 9' + n)/2 - keAR'{9'). It 
is customary to desymmetrize this map a bit by putting 
1 = 1 + kAR'{9). Then the map is 

1 = !' + keV'{9') 

9 = 9' -TT^ 2l/k (24) 

This map has been obtained by other methodsHii. 

If qualitative understanding is the main motivation, 
S may be simplified while keeping the important physics 
intact. For example, we may regard Eq.(|2^) as exactly 
defining a map which we wish to study classically and/or 
quantally. A slight but natural simplification is to use 
69 = 9' — 9. It seems natural to take AR as the sim- 
plest possible. Taking AR = cos 9 as the simplest possi- 
ble periodic perturbation ildelds the well known Chirikov- 
Taylor or standard mapcJCJ. Although this choice of AR 
does capture most of the phenomena of interest, we find 
that it is not typical in certain senses. We return to this 
point below. 

Let the invariant loop be given by a formula Iinv{9)- 
This function will have to be two valued if it is defined 
for a finite range of angle, and it can be single valued if 
it is defined over the full range of angle. It must satisfy 

T{liNv{9'),9')^ {liNv{0),9). 

In terms of the generating function, it must satisfy 
(9) = dS{9,9')/d9, {9') = ~dS{9,9')/d9'. In- 
tegrating this relation, and calling F{9) = j ljMv{9)d9 , 
we see that 

F{9) = S{9,9*)+F{9*)+const (25) 

where 9* [9) is the angle such that dS{9,9*)/d9* + 
F'{9*) = 0. 

If S{9, 9') + F {9') appears as a phase in an integral 
over 9', the stationary phase point is of course 9' = 9*. 
Thus we see that we find an approximate invariant loop, 
liNvid) « k^f'id) = ±k^y'Em-Vi9). More gener- 
ally, this is the first term in a series expansion for Ijnv- 
Therefore, Em is an approximate constant of the classical 
motion, which, like the energy, will be quantized. 

The loop liNV we have found so far is not invariant in 
a more accurate assessment. Rather, it will be mapped 
into a new loop T{lmv) h{9) ^ Iinv{(^)- The area 
enclosed by the two loops is the same, since T is area 
preserving. It is possible to estimate the phase space area 
of points inside li but outside Iinv- In the worst case, the 
area occupied by such 'turnstile' points is proportional 
to e'^/^, which is a factor e smaller than the area of Iinv 
itself. If the turnstile area is much smaller than fi, i.e. if 
kb'^ << 1, it is reasonable to believe that Iinv is a good 
approximation for the purposes of quantum mechanics. 
Of course, it may be possible to make a correction to 
liNVj i- e. to find a higher order approximation to the 
invariant loop such that the turnstile area is even smaller. 
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This would then allow us to take kb^ > 1 and still obtain 
good results. 

We thus see that the libration states with Em < 
maxl^, are to be identified with the resonance islands 
about the stable fixed points. The rotational states with 
E„i > max V are distortions of the unperturbed invariant 
loops liNV — const. If Em » max V, this distortion is 
small. However, this approximation breaks down at some 
stage because there will be higher resonances which must 
be taken into account. 



F. Energy and quasienergy quantization 

We have seen that there is a WKB quantization of the 
approximate classical constant of the motion Em- This 



leads to an expression for the quasienergy uj, of Eq.(lO) 

uj ^ Lu{k) ^ 2k + keEm +oJa +tt/2 (26) 

Table 1 gives a comparison of numerical quasienergies 
in comparison with this formula for|-a. number of states. 
The method used is that introducedEII to solve the stan- 
dard map. Starting with an approximate wavefunction, 
a long time series is obtained by applying the T opera- 
tor repeatedly. Fourier transform of this series yields the 
eigenphases and eigenfunctions. The T operator can be 
applied very efficiently if it can be factored into a part 
dependent only on 9 — 9' , and a part whose phase is a 
sum of a function of and a function of 9' . The latter 
is not strictly true for the nearly circular billiard, but 
it is true to order b^. Note also that we are comparing 
our method to the numerical solution Eq. (jsj) , and not to 
the numerical solution of the Helmholtz equation. The 
question of how well the T operatoi' corresponds to the 
billiard has already been addressecH. 

The results are quite good, that is, the errors are small 
compared with the separation of the levels. The separa- 
tion is of order 



ke{E, 



m+1 



— Ery 



(27) 



The energy levels are given by choosing values of k 
which solve 



uj{k) — 2Tm. 



(28) 



This has solutions k ~ kn.m as for the perfect circle. [For 
AR ~ 0, UJ ~ 2Tm reduces to 2k+m? / k+Trm — (n— i)27r 
where m is the angular momentum. This is equivalent 
to Debye's approximation to Bessel's function, valid for 
k large and m/k small.] 

For fixed m, the variation of uj with k from the first 
term, 2k, of Eq.(|2^) dominates. Thus we have 



kn 



(29) 



The error in the determination of a given energy level 
is not, however, small in comparison with the mean spac- 
ing in energy of all the levels. The levels so far found 
are a small fraction of the levels in a given energy range. 



There are levels belonging to larger angular momenta and 
smaller radial wavenumber is the same range. In terms 
of k, the level spacing of all the levels is about 2/fc, which 
is the order of the size of the errors committed and is the 
order of the absolute error of a given level, as is shown 
in table2. 
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1.0097 


6.1222 


6.1218 


13 


1.1223 


6.1904 


6.1908 


15 


1.1272 


6.2553 


6.2550 


17 


1.3254 


0.0317 


0.0319 


19 


1.4175 


0.0885 


0.0884 


21 


1.5040 


0.1411 


0.1413 


23 


1.5852 


0.1912 


0.1911 


25 


1.6614 


0.2376 


0.2378 


27 


1.7327 


0.2814 


0.2814 


29 


1.7989 


0.3217 


0.3220 


31 


1.8599 


0.3592 


0.3593 


33 


1.9152 


0.3924 


0.3932 


35 


1.9638 


0.4243 


0.4229 


37 


2.0010 


0.4461 


0.4458 


TABLE I. Numerical and theoretical 


quasienergies, com- 


pared for states with different m, but 


all belonging to the 


same low angular momentum resonance. 


AJ? corresponds to 


the stadium billiard. 


WKB 'energy' parameter is also given, 


k = 


1000 and e = 6.1 X 10"" are fixed. 




n 


Em + 4/7r 


k numerical 


k theoretical 


319 


1.0108 


998.3207 


998.3213 


318 


1.0128 


995.1784 


995.1789 


317 


1.0148 


992.0358 


992.0364 


316 


1.0169 


988.8934 


988.8940 


315 


1.0190 


985.7509 


985.7515 


314 


1.0210 


982.6085 


982.6090 


313 


1.0231 


979.4660 


979.4666 


312 


1.0252 


976.3235 


976.3241 


311 


1.0273 


973.1810 


973.1816 


310 


1.0294 


970.0387 


970.0392 



TABLE IL Energies k with different quantum numbers n, 
but same m = 11, computed numerically and found solving 
Eqs. (||) and (|2| ). e as in Table 1. 



It would be nice to have a theory giving the energy 
levels with absolute accuracy less than thej-mean level 
spacing, of course. Primack and SmilanskyEa have dis- 
cussed how this might be possible within the framework 
of quasiclassics, even though the errors made in arriving 
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at the T operator and doing stationary phase integrals 
are of order Ti. This however involves finding the mean 
density of levels to better than leading order accuracy by 
a separate calculation, and combining it in a particular 
way with sums over periodic orbits. 

Nevertheless, the results we have obtained are very 
useful and contain nearly all that is desired. The en- 
ergy levels are classified into groups. Matrix elements 
of smooth operators are large only between levels in the 
same group. The error is small compared with the spac- 
ing of the levels in the same group. The overall statistics 
of the levels, when found on the scale of the mean level 
spacing, is also given correctly. Long rpige correlations 
of energy levels are also given correctlytll. 



G. Wavefunctions and localization 

The 'wavefunctions' are also given to good approxima- 
tion, as we showed above. 

Of course, there could be an accidental degeneracy be- 
tween energy levels coming from small angular momenta, 
such as we have discussed above, and energy levels car- 
rying large angular momenta, which we will calculate 
below. Then, the true eigenstates will be some appre- 
ciable admixture of large and small angular momenta. 
Lazutkin's rigorous results also allow for this possibilityQ. 
However, the matrix elements of the T operator between 
widely different angular momentum states is exponen- 
tially small. This means that such accidental degenera- 
cies will be rare. 

In the presence of symmetries, there will be degen- 
eracies which are exact in the absence of the exponen- 
tially small coupling. An example is the states asso- 
ciated with the periodic wells discussed above. In the 
presence of time reversal invariance, the states associ- 
ated with positive angular momenta are degenerate with 
those of the corresponding negative angular momenta, 
again with very small splitting due to tunnelling between 
the two momentum regions. If desired, this splitting can 
be estimated within the framework of oup-theory. This is 
not entirely trivial as it requires a studytil of 'resonance 
assisted tupnelling' which is an analog of 'chaos assisted 
tunnelling O. 

The wave function so far found contains only small 
angular momentum components. This can be seen ex- 
plicitly in angular momentum representation, determined 
schematically from the integral 

iJi= J dO {Em - V{e))-'^ exp {ik^ef{9) - iW) . (30) 

If /jy^ >> 1, stationary phase may be employed. Then, 
it is found that only the neighborhood of the angles solv- 
ing 

k^ef'id) = kVeVEm - V{9) = I (31) 

contribute appreciably. A range of I centered roughly at 
k^/ eEjn and of width of order fc\/e maxV can satisfy this 



condition, so fcV e max V is a sort of localization length 
for this kind of wave function. Note that only librational 
states of the lowest resonance have much overlap with 
zero angular momentum. 

Outside this range, the angular momentum com- 
ponents of the wave function decay exponentially for 
smooth V. For the stadium case, V = |sin6'|, it may 
easily been seen that the decay is as rather than ex- 
ponential. This makes it necessary .tci, iLse. nonstandard 
definitions of the localization lengthE30'L3 if statistical 
results are to be calculated. 

Carrying out the integral of Eq.(|30|) gives an expres- 
sion 

^l^Y. I^'(^a)r^ expz (k^efiOa) - Wa + Va\) (32) 

So 

where 9a are the solutions of Eq. (|3l| ) . 

Notice that, if as expected, V has different signs 
at the different 6*^, there will be an additional Maslov 
phase index Va distinguishing the two solutions. Also 
notice that the singularity of the prefactor at the classi- 
cal turning point in angle, where V {9) = Em has disap- 
peared. Instead, the angular momentum regions which 
make V' vanish, which are the turning points in the an- 
gular momentum representation, are singular, and should 
be treated by a technique going beyond the first WKB 
approximation. This shift of the region of breakdown of 
the leading WKB approximation under Fourier transform 
is the basis of Maslov's treatmentcHi of this subject. 

The solutions thus found are localized in angular mo- 
mentum space. From the point of view of KAM theory 
and perturbation theory, this is unremarkable. First, the 
KAM theory predicts that the classical motion will be 
confined to much the same region. Second, matrix ele- 
ments of the T operator in the angular momentum rep- 
resentation are small away from the diagonal. In fact, 
consider Tw = (27r)-i/ J d9d9'e-'^'^T{9,9')e'^'^' where 
{I — l')/k is of order unity. The 9' integral will be sta- 
tionary near dS[){9 — 9')/d9 = I', and the 9 integral will 
be stationary near dSo{9 — 9')/d9 — I. Restricting con- 
sideration to the standard case where there is only one 
stationary point, we see that the integral will be expo- 
nentially small. 

Such a band jdiagonal matrix Tw is studied in lo- 
calizatpj^-jjliypi-yEZl and is effectively the subject of the 
papersOEBBj^EJ cited earlier. Of course, since the long 
time behavior is diffusive, in the non KAM cases, it was 
not obvious exactly the mechanism for localization. If 
stationary phase arguments are sufficient to get the width 
of the band, as we have assumed, the problem is solved. 
The cited papers are concerned with the non KAM case 
that ke^ is sufficiently large that stationary phase breaks 
down. 

It should be emphasized that our whole theory is based 
on the assumption, validated self consistently, that the 
eigenstates are localized in angular momentum. 

Although it is usually not worth the effort, the full 
two dimensional wavefunction can also be found in WKB 
approximation. 
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H. Large angular momenta 

We turn to a discussion of the problem of larger angu- 
lar momenta. As we have just remarked, a large angular 
momentum state cannot communicate directly via the T 
operator with its negative angular momentum counter- 
part. It's true that time reversal symmetric systems will 
be symmetric under / —I, so that true eigenstates will 
be a mixture of these states. However, the even and odd 
eigenstates will be practically degenerate, and we can 
consider one sign of angular momentum only to get the 
main effect. 

Suppose we wish to find eigenstates centered near 
angular momentum I. One way is to reduce the prob- 
lem to the one just solved by making a sort of gauge 
transformation, in which the eigenstates -0 e^''^tp, 
T ^ e^^'C^^'^'^f', and ftp = 'ip. We may then assume 
that "0 is slowly varying. 

There are some complications in this case as compared 
with the case where / = 0. First, to leading order, the sta- 
tionary phase point is given by dSQ{9—9')/d6 = I, namely 
at 6 — 9' = —8;, say. However, 8; has no particular rela- 
tionship to 27r, that is, it will be an irrational multiple of 
2tt, unless I is specially chosen. (The negative sign in the 
definition of 8; is purely for orthographic convenience.) 
If / is an integer, this cannot happen except at special 
symmetry values such as Z = 0. 

There are two possibilities. Either 8; is very close to a 
rational multiple of 27r, Qpq — 2'Kp/q, where q is not too 
big, or it is sufficiently far from such a multiple. It is the 
job of the theory to tell us what 'very close', 'not too big' 
and 'sufficiently far' mean in detail. The resonant case 
is the one where 8/ is 'close' to a rational, the KAM-set 
case has 8/ 'far' from a rational. 



1. Resonant case 

The first case we deal with is that of a resonance of 
period q, which involves destruction of the original in- 
variant tori and their replacement by new tori with a 
different topology. These are represented on the surface 
of section by a period q island chain. The second, nonres- 
onant, case deals with invariant tori which are not much 
distorted by the perturbation. 

In this case we consider states whose angular momen- 
tum components are not too far from Ipq where the lead- 
ing stationary phase point corresponding to Ipq is Qpq. 
Note that Ipq is not generally integer, so that ip will not 
be 27r periodic. 

The Ansatz made before. 



= exp{ikbf{e)) 



(33) 



does not lead to a solution. Indeed, expanding T, / about 
8/, we do the 9' integral and find 



exp 



ike (cf{9 + Qpq)^ + diAR{9) + AR{9 + Opq)) 



e-'"e'" exp(zfc6(/(0) - /(0 + 8^,)). 



(34) 



Here c is determined from the second deriva tive of Sg . 
The coefficient d for the billiard is d = = a/I — P/k^. 
a is the phase obtained by doing the 9' integral. For the 
case defined by Eq.(|), a = ^n + 2y/k^ - 1^-21 cos'^ij:). 
In this case Qpq = —2 cos~^ (l/k), for positive I. To satisfy 
Eq.(^^ we must require 



fi9 + Qpq) = fi9) + const 



(35) 



and 



c^if'f + diAR{9) + AR{9 + Qpq)) ^ const. (36) 

These two conditions cannot in general be simultaneously 
fulfilled. Eq.(^ implies that /' is g-periodic, i.e. peri- 
odic with period Qpq, but only for special choice of AR 
will the terms in AR in Eq.(^) be g-periodic, if g ^ 2. 
We therefore improve the Ansatz and take 



i>^exp[ikd{bh{9) + b^f2{0) 



(37) 



Introducing d, characterizing SQ{Qpq), simplifies later 
formulas. Treating 6^/2 as slowing varying, in this or- 
der of the calculation, we find the conditions 



f i{9 + Qpq) ^ fi{9) + const 



(38) 



and 

f2i9) -f2{9 + Qpq) = /f + AR{9) + AR{9 + Qpq) - 

(39) 

where Em is a constant to be determined. 

Let us define the q-average of a function of angle by 
Pqi^) = I ELi F{9 + rQpq). The Fourier components of 
Fq are just those of F with all components whose index 
is not divisible by q discarded. Making a q-average of 
Eq.(|9|) leads to 



f[^±jE,n-Vq{9) 



(40) 



where Vq = 2ARq. 

The function /2 cannot be completely determined at 
this level, since an arbitrary q-periodic function can 
be added to it and not affect Eq.(|3^). Calling Vpq — 
AR{9) + AR{9 + Qpq) ~ Vq we can solve for the re- 
maining part of /2. This function has no Fourier com- 
ponents whose index is divisible by q. The solution of 
f2{9) — 12(9 + Qpq) = Vpq can be made by Fourier trans- 
form of the equation. An explicit solution can also be 
written 



f2iO) 



1 



' + rQpq) 



(41) 



r=0 



This suggests making the following transformation. 



Let -0 



Jkudf; 



^"0 and consider the corresponding oper- 



ator T which may be approximated by 



f{9,9') 



.[Ile^p{2tkd-—S9^) 
V 27ri ^ 4 ' 

X exp{ikde{AR{9) + AR{9 + Qpq) 

+ f2{0 + Qpq)-f2{0))) 



(42) 
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The second exponential is exp{ikdeVq). If kde is small, 
we expand as in Eq.(|3|) above. Thus, tp solves the 
periodic potential problem with potential Vq and effective 
Planck's constant Xjkd^. 




0.937 



FIG. 5. Approximate, to order e, and exact invariant loops 
at the fourth order (b) and twenty-eighth order (c) resonances. 
The value of e is 0.05 in (a), (b) and 0.01 in (c), not very 
small. T] — 0.7 in (a), (b) and 0.19 in (c). Thick lines are ex- 
act invariant curves, thin lines are theory. In the fourth order 
case we have verified that the failure of the approximation is 
largely due to the absence of the order fo^ contributions - in 
(a) theoretical curve including b"^ order is plotted. Higher pe- 
riod resonances have relatively, although not absolutely, big- 
ger corrections from higher orders in y^. 

If kde is of order unity or larger, the argument goes 
through as before, so that we have obtained the solution 
to the resonant case. This also gives the pref actors, as 
before. 

We show in Figs. 5a, b 5c some invariant loops for 
higher resonances pq = 1,4 and pq — 3,28 respectively. 
The numerical loops obtained by iterating the SS map are 
compared with the theoretical Ijnv — Ipq + bf[ + 6^/2 • 

There are several things of note. First, the contribu- 
tion of /2 becomes relatively more important for higher q. 
This is for two reasons. Most obviously, (eVq)^^'^ system- 
atically becomes smaller with increasing g, and eventu- 
ally becomes smaller than e, so that the small resonance 
islands 'float' on a relatively large 'wave' supplied by /2. 
Second, /2 can itself become rather large. This can be 
seen from formula (41), which can become fairly big if 
the terms in the sum add 'in phase'. This effect is most 
pronounced near a large resonance island, e.g. at g = 2. 
Then the term e/2 for p/q slightly greater than ^ with q 
large, must be almost as big as \/ef[ for q = 2. 

Another remark is that higher order terms more gen- 
erally are relatively more important for large q. This can 
be seen in the results (Fig. 5c). These deviations are 
apparent in this case of rather large b = 0.1, 6^ = 0.01, 

= 0.001. Again this is because Vq is becoming smaller, 
while the higher order corrections presumably remain 
roughly constant. The observed deviations are appar- 



ently of order . Of course, if kb^ is small, the quantum 
results won't be affected by the deviations. 



0.2 
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FIG. 6. The periodic quantum state for the resonance 
pq — 1,3. Theoretical wavefunction is cos{kdbfi) ■ 
sin (^W + kdb'^ f2)\ in (a), (b) we integrated out exact and semi- 
classical functions locally with respectively sm (16), cos [16). 
k = 4032, e = 6.79 x 10~*. 

The case pq — 3, 28, shows an impressive set of wrig- 
gles coming from /2 which are quite well predicted the- 
oretically. We have not tried to find the 'best' value Em 
corresponding to a given classical numerical Iinv, ex- 
cept by trial and error. This may contribute somewhat 
to the discrepancies. The islands are of different shape, 
however, and vary in height from one island to the next. 
This effect cannot be reproduced to the order which we 
have worked. 

In Fig. 6 we also show the quantum states correspond- 
ing to the resonance = 1, 3. What is plotted is 1/;, that 
is, the rapidly varying factor exp il9 is removed. 



2. Quantization conditions 

The quantization conditions are somewhat modified, 
since the central angular momentum Ipq is not usually 
integer. We require that exp{ilpq9 + ikdbfi[0)) be 27r- 
periodic. For the rotational levels, E,, 
a modified quantization condition 



> max Vq we have 



/ de^/E~V = 27r(r 
Jo 

where —5 is the fractional part of Ipq. 
If there are solutions of 



5)/kbd 



(43) 



kbd 



deJEm~Vq{e) =n{m + i^) 



(44) 



where the integral is between turning points straddling 
a minimum of the potential, v again being the Maslov 
index, there will be librational solutions repeated q times. 
There are then q nearly degenerate states, and the phase 



13 



shift from one well to the next will be slightly different 
than Qpq in such a way as to make the total wave function 
27r-periodic. 

These quantization conditions make the area of the 
loops an integer number of Planck's plus the correction 
due to Maslov indices. Thus, it is equivalent to EKB 
quantization. 



3. When the resonant solution is applicable and needed 

The crucial scale is usually provided by Vq . We assume 
without loss of generality that AR and thus Vq have van- 
ishing mean. 

If librational states exist, then we clearly need the res- 
onant solution. This requires Vq > q^/k'^e. 

If Vq is very small, the smallest rotational Em is of 
order (kbd)'"^. Thus, if Vq > {kbd)~^ we cannot 

neglect it. If Vq is much less than this value we can ne- 
glect it and the states in this angular momentum neigh- 
borhood are candidates to be treated by methods valid 
for the KAM set. 

We remark on several cases which are distinguished by 
the way the Fourier components of AR behave. 

If AR is a pure low harmonic, or sufficiently close to 
one, e.g. AR = cos 9 or AR — cos 20, then Vq = for 
(7 > 2. Of course, it is immediately recognized that these 
cases are to first approximation a shifted circle and an 
integrable ellipse, respectively, so that they can be trans- 
formed to a perturbation about an integrable system with 
perturbation parameter rather than e. 

From our perspective, if Vq vanishes or is sufRciently 
unusually small, one should determine if some higher 
power of e gives resonant behavior. This is not hard 
to do in the simplest cases, but is tedious for more com- 
plicated cases. For example, the order q resonance of 
the standard map for small e will have width parameter 
proportional to e''/^ , since the effective resonant coupling 
is (ecos6')^. Cases where the Fourier expansion of AR 
is a finite Fourier series will have a rather complicated 
condition which is difficult to treat in generality. 

Analytic perturbations usually have Fourier coeffi- 
cients which drop off exponentially with index. This leads 
to an estimate, 77^. If 77 > e, we can expect that the 
treatment given above is adequate and the leading or- 
der expression for the resonance is proportional to -^/e. If 
not, we either have to work out the higher order effects or 
take k sufficiently small that we can ignore the resonance 
entirely. 

Finally, there are the cases with nonanalytic AR, 
whose Fourier coefficients drop off as a power, say. The 
billiard AR = \smd\ has Vq ~ q~^, for example, where 
KAM theory breaks down completely. We can generally 
expect our leading order solution to be valid for these 
cases. In this case, the entirep^hase space is best re- 
garded as filled with resonancesEj, leaving no KAM set. 



4. Close to a large resonance 



The following issue is to some extent still unresolved. 
Suppose Vq is sufRciently large that Eq.(|4^) holds non- 
trivially. Consider, for example, the resonance labelled 
pq — 7,16 in Fig. 2. For values of Em > maxVq 
the predicted invariant loops are well described by the 
low order resonance formula, although they don't give 
higher order resonances. Assuming Em large enough to 
expand in Vq, the invariant loop is given by bf{ + b^f2 ~ 
by/Emi^ ^ \Vq/Em) + &^/2- This may be shown to ap- 
proximate the nonresonant solution, also called /2 but 
for a nonresonant value of I, found in the next sectionEil. 
However, even if one is quite close to the separatrix, 
where one expects the resonant formula to be best, the 
resonant or nonresonant formula at the appropriate I 
value is not bad. Even the formula for a higher order 
resonance which seemingly lies in a region strongly af- 
fected by the big resonance does qualitatively quite well, 
although there are some small discrepancies. 

This seems to present a contradiction. Close to the 
big resonance, e.g. p — \,q = 2, a, secondary island chain 
from a high order resonance follows well a curve given 
by b\/ Em — V2 and also a curve with seeming smaller 
variation e/2, corresponding to a resonance pg. However, 
according to Eq.pl|), /2 can become quite large for large 
q, which can compensate for the fact that e << ^/e. 



5. Nonresonant case, quantization of the KAM set 



In this case, we start with a value of integer angular 
momentum I, such that O; is not resonant. There is no 
term /i and the solution begins with /2. The condition 
for /2 is 

f2{0 + Qi) - hiO) = ^R{0 + Qi) + AR{9). (45) 

This may be solved in terms of Fourier components of 
AR , denoted AR r. The result is 

1 fD^rQi I 1 

M')-^J:-I^,^^re-'- (46) 



This sum will be well behaved, provided ARr, drops 
off sufficiently rapidly, and is small when the denomi- 
nator is small. The denominator is small, when r = q 
and 0; is close to a rational Qpq, with q not too large. 
The denominator then approximates i{qOi — 2'iTp). If this 
condition does not hold, the resonance treatment of the 
previous sections must be employed. 

At this order of the calculation, there is no shift in the 
energy levels for the KAM set states. 
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FIG. 7. The golden mean invariant torus, approximated to 
order e [solid line cos ^ +6/2(0)], and numerically exact [clas- 
sical map represented by dots; 20,000 mappings]. Parameters 
are e = 1 x 10"*, i] = 2 x 10~^. The inset enlarges a cusp- like 
portion of the figure to display the degree of failure of the 
approximation. This torus has Oi = 27rw — 7r(\/5 — 1), and 
it is 'farthest' from a low order rational. Such a torus is not 
really smooth, and is not describable by a convergent power 
series. It is 'transversally smooth', which is enough to control 
mathematically its propertiesu. One can pick out the closest 
resonances, e.g. 5/8, 8/13, etc. given by the Fibonacci series. 
Presumably, as rj is decreased, this is the last torus remaining 
in the KAM set. Note that the new tori at the resonances are 
relatively smooth. Presumably, on the fine scale of e^^' they 
would also have such a structure. 



quantum state that corresponds to the same KAM set 
was build and compared to the numerical wavefunction 
in Fig. 8. 



I. Higher order calculations 

We briefly mention how one goes to higher order in 
the calculation. Some details are put in the Appendix. 
We confine attention to the resonant case. 

Basically, one makes the Ansatz 



V' = exp [i (^Ip 



k{bh 



•))] (47) 



Instead of expanding Sq to second order m 69 ^ 9' — 6 — 
^pqi fi to first order, /2 to zeroth order, and neglecting 
/a, we expand, step by step, Sq to the A/'th order, /i 
to the M — I'th order, /m to the zeroth order, and 
neglect /m+i, where M > 2, and kb^^~^^ << 1. The per- 
turbation, which we may call eS2{9, 9') + t^Si{9, 9') + ... 
is expanded according to the same rule. The order r cal- 
culation determines the non-periodic part of fr and the 
periodic part of fr-i- 

The resulting integral is not done exactly, but by sta- 
tionary phase in the neighborhood of 59 — 0. The sta- 
tionary point is regarded as shifted by AOp^ which can 
be expressed as a power series in b. The b^ term in this 
shift is bigger than the width, l/Zc, of the effective region 
in the 59 integral if >> 1. A sample calculation is 
done in the Appendix. 

There are interesting questions about how this method 
breaks down if there are unusually large derivatives of 
Ai?, but we shall not address them in this paper. 
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III. OTHER INTEGRABLE SYSTEMS 
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FIG. 8. Exact and approximate wavefunctions for the pa- 
rameters of Fig. 7. The factor e*'" has been removed and the 
real part of the wavefunction is shown. 

Formula (4^) has been tested numerically. In Fig. 7 we 
plot /2(^), the derivative of Eq.(|4^) and compare it with 
the numerical invariant loop, obtained by propagation 
of an orbit. It is apparent from this that the invariant 
loops are not totally smooth and featureless. KAM the- 
ory controls the singularity of these curves, although they 
cannot be represented by a convergent power series. The 



There are an infinite number of integrable systems 
whose perturbations can be studied. Another interest- 
ing billiard is the rectangular billiard. As an example of 
a perturbed rectangular billiard, consider the trapezoidal 
billiard, whose sides consist of the a;— axis, < 2: < a, the 
vertical segments, 0<y<6 — e, {x = Q) , Q < y < b + 
[x = a), and the slightly sloping segment connecting the 
upper ends of the vertical segments. 

The results arc much the same as before. We just 
mention what is found for tVq. This potential is periodic 
in X, with period 2a /q. In one period, the potential is an 
isosceles triangle 'tent', of zero mean, and peak absolute 
value e/<7^. Note that x is continued outside the physical 
region. 

In fact, for any reflection symmetric perturbation Ai? 
of the circular billiard, there is a corresponding perturba- 
tion of the rectangular billiard, in which only one side of 
the billiard is perturbed. The leading orders are almost 
the same in the two problems, except that the param- 
eter d, having to do with the second derivative of the 
integrable action Sq is unity in the rectangular billiard. 
Higher orders are different, as well, which may be helpful 
in studying those effects. 
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IV. NO SMALL PARAMETER 

In this section we briefly take a case where some 
states of a system can be found, even though there is 
no small parameter. Consider the standard stadium bil- 
liard, which has a straight side of length 2a, as in Fig. 
9. The side a is arbitrary, but not too small. We want 
to study the 'bouncing ball' states which have low linear 
momentum parallel to these sides, but a high perpendic- 
ular momentum. 



Em is sufflciently small, then / varies slowly, and the ex- 
pansion is justified. The possibility of an emergent small 
parameter, such as Em allows us to use a perturbation 
approach. 

It's clear that a sufficient number of classical orbits 
giving a localized wavefunction is a precondition for spe- 
cial solutions of this kind. This can happen if the classical 
orbits are nearly nonisolated. That will happen system- 
atically near stable periodic orbits. 

Some further remarks are made in an Appendix. 
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FIG. 9. Half of stadium billiard with long straight side. 
Notations as in the section iv. 

Reduce the billiard to its upper half by symmetry. 
Take as nominal integrable system the infinite channel 
of width i?o = 1 as shown. We take as SS the upper 
boundary of the billiard, and use position x, measured 
from the symmetry point, on the upper channel to label 
position on the surface of section. The billiard boundary 
is described by distance ^ from the upper nominal chan- 
nel, i{x) = 0, \x\ < a, ^{x) « (|x| - a)V2, > a. The 
integrable action is, for states odd under reflection about 
the horizontal symmetry line, 

^"^"""^'^ = 2kJl + hx-x')^ ^2k+%- x'f. 

(48) 

States even about this line have the same action with 
an additional term nh, that is, the T operator has an 
additional overall negative sign for the even states. The 
change of sign of the second derivative, as compared with 
the circle, means that the attractive and repulsive regions 
of the potential are reversed in comparison with the pre- 
vious case. We can approximate S2 by 

S2ix,x')^-k{a^)+a=c')). (49) 

Put 

V{x) = S2{x,x) = 0; \x\ < a 

= {\x\-af]\x\> a (50) 

We look for a solution 

V' — exp{ikf{x)) (51) 

where / is slowly varying. One finds, upon approximat- 
ingf{x')^f{x) + nx)ix'-x) 

fix) - ±^E~Vix) (52) 

Thus, we look for a solution of a particle in the potential 
well V{x), and with effective Planck's constant 1/k. If 



V. SUMMARY 

Weakly perturbed integral Hamiltonian systems in two 
dimensions have been much studied. Up to now, the 
methods employed, although effective numerically, did 
not lend themselves to visualization and simple formulas. 
By using surface of section techniques as generalized to 
quantum systems by Bogomolny, we obtain results very 
much like familiar results from WKB theory, at least in 
the leading orders. 

The relation of the perturbation parameter e, to the 
dimensionless wavenumber k, or equivalently to 1/h is 
important. The size of fee, the number of wavelengths in 
the change of a billiard boundary, is not especially impor- 
tant, for resonant states. Rather, the main parameter is 
k^/e , which determines the Shuryak border demarcating 
the parameter range where ordinary quantum perturba- 
tion theory applies and where it fails. Results are quite 
simple if ke^^^ is small. In that case, the derivatives of 
the perturbation do not enter directly into the solution. 
For nonresonant states, the size of ke determines whether 
ordinary perturbation is valid or not. 

The results can also be regarded as giving expressions 
for the invariant loops, the intersections of invariant tori 
with the surface of section. Knowledge of these, rather 
than of the whole invariant torus structure is sufficient 
to quantize the system, quasiclassically. 

VI. FUTURE WORK 

Several avenues for future work are underway. It would 
be interesting to see how the perturbation approach 
breaks down in some detail, in cases such as the stadium 
billiard. Working out and verifying higher order effects in 
detail will be of some interest. It would be well to extend 
the method to arbitrary surfaces of section rather than 
the special ones we have used. A number of interesting 
problems are raised by cases such as AR — cos 29. Other 
perturbation systems would be intepjsting to study, for 
example, coupled quartic oscillatorsE3. It would be well 
to extend the method to weakly coupled harmonic oscil- 
lators. A general theory of whispeciag gallery orbits in 
sufficiently smooth convex billiarda^Ej, might be possi- 
ble, but would require a better T operator than we have 
used above, since diffraction effects become important 
for short classical skips between nearby boundary points. 
Extension to d > 2 may be in principal possible but in 
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practice is likely to be difficult. There are interesting 
numerical results, for example in ray splitting billiardsE3 
that suggest it might be possible to solve the case that 
two or more orbits contribute to a specific T operator, 
requiring a sort of matrix version of our theory. Finally, 
it could be hoped to use this technique to studjy some 
or even most states in pseudointegrable systemsEll, which 
have nonisolated orbits which however do not lie on a 
simple toujs in phase space. Many of these projects are 
underwaytil. 
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APPENDIX A: 1 < Ke^^'^Ke^ « 1 

We here indicate how to find the resonant solution if 
kb^ ~ l,fc6^ >> 1, and kb^ « 1. We expand So about 
the nominal stationary point Qpq to third order: 



The last four terms are of order 6^, while the first two are 
of order b, since the leading expression for 69 is of order 
b. Thus, we find 

60 = 2bf[ + 2b\2Adi{f[f + 2f[f'{ + + Vi) 

The integration is performed, yielding a condition at or- 
der b^ 



di/f + 2/f/r + 2/{(/^ + y2') 
= h{0) - W + Qpq) + c^onst 



(A6) 



^o( 



9') + lpq{e' -e)~ SoiOpq) + IpqQp 



-k-6e^ + lkddi6e^. (Al) 
4 8 



Taking the q-average of both sides, we have an equation 
determining /2, the part of /2 invariant under q-average. 

f2, = -K~fifi-ldi{.f[r + jj (A7) 

The part of /a vanishing under g-average is determined 
from 

f3{0) -h{0 + Op,) = 2f[(f!, - + ~ K,',) 

and /a is determined by the next order condition. Note 
that if /i is double valued, so if f^. It can therefore change 
the shape of the resonance islands. 



APPENDIX B: NO SMALL PARAMETER 

We add some further remarks concerning the quanti- 
zation of the bouncing ball states in the stadium billiard 
and some related systems. 

A WKB quantization condition of Eq. (^2|) is 



(A2) 



2k 



I dx^/Wry 

Jo 



Vix) 



The order e term is 

52(6*, 6*') ~ kb^{dV2{e) + dV^{9)6e) 
and the order term is 

S4{0,e') = kdb'^V4{e). (A3) 

We have normalized to the second derivative of Sq which 
is written in terms of the constant d. Note that V2 is not 
the derivative of V2 . 
The wavefunction is 

^ = exp(-5) expikd + 6^/2 + b^fs + ^^A)) • (A4) 

We expand 

That is, he takes 1^ = and neglects the integral for 
fi[9') ~ fi{9 + Opq) + 69f[{9 + Qpq) + 69^ fi{9 + Opq) \x\ > a. He considers just the states odd under refiection 

X — !■ —X, which for us means m odd, i.e. m 2m — 1, 
E,n = (7r(m- l/4)/fcaf . 

The equation for the energy of the odd-odd states is 
then 



The term in Em comes from the integral in the region 
X > a andjs smaller than the leading term. 

TanneitJ, by a method using only periodic orbits and 
with a different surface of section than we use, finds after 
an impressive calculation, an approximate energy equa- 
tion equivalent to 

r— 7r(m + l/2) 



2ka 



f2{9') = f2{9 + Qpq) + 69f^{9 + Qpq) 
h{e') = h{9 + Qpq) 



(A5) 



We know /i and the part of /2 which has no Fourier 
components which are multiples of q. We investigate the 
conditions on the solution which arise at order kb^ . 

The stationary phase point of the integral satisfies the 
equation 



1 



bf[ + 2Ui69'^ + bf'{69 + 9 f!^ + b^V^ = 0. 



Tanner finds this to be a good approximation, although 
he applies it outside the range of our derivation, that is, 
for a = 5, he takes n — 1, 2, 3 and m = 1..9. 
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We find things a bit more complicated, especially in 
determining the Maslov exponent. This will be reported 
separately, however. 




FIG. 10. The ice cream cone billiard with parameter 
/3 = 7r/4. 

Next, consider an ice cream cone billiard, which is a 
unit circle for \6\ > (3, and a triangular shaped region for 
\0\ < (3. See Fig. 10. Microwave expempents on such a 
billiard have been recently carried outO. The equation 
for the straight line segments is 

\y\ — CSC (3 ~ X cot (3 

for cos (3 < X < sec (3. 

Using the unit circle as nominal integrable system, the 
billiard boundary as SS, and as the parameter deter- 
mining position on the SS, we find for the action, on and 
near the circular section. 



S{e,e') = 2fcsini(6l- 



Km+m) 



where ^{9) » i(|6i| - (3f- \e\ < /3, = 0; /3 < |6i| < tt. 

This leads to essentially the same WKB problem as in 
the stadium case just discussed. The change of sign of the 
second derivative is compensated by the change of sign 
of the deviation from the integrable system. Because, 
however, zero angular momentum makes 0' « ± tt, the 
effective potential V2 = \ [\e\ - (3f; \9\ < (3, i(|7r ^ 9] - 
\tt - 9\ < (3, = 0; P < \9\ < TT - (3. In other words, 
there are two 'potential wells' whose confining walls rise 
at half the rate of the stadium case. 

The potential formulation seems to hold in the sta- 
dium case, even for a —> 0, or in the ice cream cone case, 
for (3 ^ TT. In fact, in the latter case, the ice cream cone 
billiard approaches a very long stadium. It is a good 
exercise, however, to show that in these cases, the solu- 
tion is such that the implicit small parameter, Em can no 
longer be taken as small enough to justify the procedure. 

Of course, it is not difficult to extend these results to 
deal with cases having an additional small perturbation. 
For example, the flat sides of the stadium might be given 
some small wriggles. An example of this type that has 
been published^ makes the radius at one end of the sta- 
dium billiard slightly larger than at the other end, so that 
the stadium flat sides are not quite parallel. 
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